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Abstract 



^ ■ Various classical counterparts for the two-level pairing model in a many-fermion 

system are presented in the Schwinger boson representation. It is shown that one of 
the key ingredients giving the classical descriptions for quantal system is the use of 
the various trial states besides the su(2) Cg> sw(2)-coherent state, which may be natural 
selection for the two-level pairing model governed by the su(2) ® sw(2)-algebra. It is 
pointed out that the fictitious behavior like the sharp phase transition can be avoided 
by using the other states such as the sm(2)®sm(1, 1)- and the su(l, l)(g>su(l, Incoherent 
states, while the sharp phase transition appears in the usual Hartree-Fock-Bogoliubov 
and the quasi-particle random phase approximations in the original fermion system. 
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§1. Introduction 



The use of boson representation for many-fermion systems gives a powerful method to 
describe the dynamics of the original fermion systems. On the basis of the Lie algebraic struc- 
ture of the model Hamiltonian, the Schwinger 1 ) and the Holstein-Primakoff 2 ^ type boson re- 
alizations give helpful tools to investigate the dynamics of the original fermion systems. The 
boson mapping method, which was proposed in Refs.3), 4) and 5), is also useful method 6 ^ to 
describe the many-body systems beyond the usual mean field, the Hartree-Fock(-Bogoliubov) 
and the random phase approximation. The Marumori-Yamamura-Tokunaga (MYT) boson 
mapping method 4 ^ is characterized by the mapping of operators acting on the Fock space 
in the original systems into operators acting on the boson Fock space in the corresponding 
boson systems. This idea was strongly stressed by the present authors and widely applied 
to the many-fermion and many-boson systems to study the dynamics under consideration. 7 ^ 
The oscillation around the mean field configuration is, of course, described in terms of the 
bosonic degrees of freedom. 

As another advantage of the use of the boson representation of the original fermion 
system, we can say that it is possible to introduce various trial states in the framework 
of the time-dependent variational method. In the fermion system, the trial state may be 
restricted to a simple class in the mean field approximation. For example, if the fermion 
system is governed by the su(2) <E> sw(2)-algebra constructed from the bilinear forms of the 
fermion creation and annihilation operators, the trial state in the variational method may 
be restricted to the Slater determinantal or BCS state, namely, the su(2) ® sw(2)-coherent 
state in the framework of the mean field approximation. This is just the case in which the 
two-level pairing model in a many-fermion system is considered. On the other hand, in the 
boson representation of the original fermion system, various trial states are possible, even 
if the system is governed by, for example, the su(2) <g> sw(2)-algebra. The su(2) <g> su(2)- 
coherent state corresponds to the Glauber coherent state in the boson representation. In 
addition to this state, the su(2) <g> su(l, 1)- and the su(l, 1) (g> su(l, Incoherent states are 
possible to describe the dynamics of the original two-level pairing model governed by the 
su(2) <S> sw(2)-algebra. Thus, it is possible to give the various classical descriptions for the 
original unique many-fermion system. 

In this paper, we deal with the two-level pairing model governed by the su(2) <g> su(2)- 
algebra. This model is treated in the Schwinger boson representation for the su(2) ® su(2)- 
algebra in which four kinds of boson operators are introduced. In order to give a classical 
description of the pairing model, three types of trial states are introduced, namely, the su{2)® 
su{2)-, su(2)®su(l, 1)- and su{\, \)®su{l, Incoherent states. The su(2)®su(l, Incoherent 



2 



state for the two- level pairing model has been defined in Ref. 8), which is referred to as (I). 
Also, the su(l, l)®stt(l, Incoherent state for the two-level pairing model has been defined in 
Ref. 9). In this paper, we numerically calculate the ground state energy and the frequency of 
the small amplitude oscillation around the energy minimum state, and compare them with 
the exact results. Then, it is pointed out that the fictitious behavior like the sharp phase 
transition can be avoided by using the su(2) ® su(l, 1)- and the su(l, 1) <8> su(l, Incoherent 
states, while the sharp phase transition appears by using of the su(2) <g> sw(2)-coherent 
state which corresponds to the usual Hartree-Fock-Bogoliubov and the quasi-particle random 
phase approximations in the original fermion system. As an additional remark, it should be 
mentioned that, from a viewpoint different from the present one, the authors have already 
reported numerical results for the two-level pairing model, which supplement the present 
results. 10 ^ 

This paper is organized as follows. The outline of the two-level pairing model and its 
Schwinger boson representation with the four kinds of boson operators is presented in the 
next section. In §3, re- formation of this model in terms of one kind of boson operator and 
its classical counterpart is given and three types of boson coherent states are introduced. A 
classical description based on each one of the three types of coherent state is given in §4. 
In order to investigate the ground state energy with quantum correction and the frequency 
of the small amplitude oscillation around the static configuration, the quantal treatment is 
given in §5. The numerical results are also shown in this section. The last section is devoted 
to concluding remarks. In appendix A, some results in the case of the su(2) ® sw(2)-coherent 
state are supplemented. In appendix B, the derivation of quantal Hamiltonian is given in 
order to introduce the ground state energy with quantum correction and the frequency of 
small amplitude oscillation around the ground state. 

§2. Outline of the model 

In this section and partly in the next one, we list some of the relations appearing in (I). 
They help us to understand various results shown in this paper. The Hamiltonian of the 
system which we intend to describe is given in the relation (I-3-3a): 

H = e (S„(+) - So(-)) ~ G (X(+) + £+(-)) (S-(+) + S-(-j) ■ (2-1) 

Here, e and G denote the energy difference between the upper (specified by a = +) and the 
lower (specified by a = — ) level and the strength of the interaction, respectively. The set 
(S±fi(cr); a — ±) obeys the s«(2)-algebra and it can be expressed in the form shown in the 
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relation (1-3-1): 



S+(a) = M*X , S4<t) = hb* a a a , S (cr) = {h/2){a* a a a - b%) . (2-2) 

The operators (a CT ,a*) and denote boson annihilation and creation operators. The 

form f!2-2p is identical to the Schwinger boson representation of the sw(2)-algebra. From the 
composition of the Hamiltonian (12-11) . we can understand that the present system obeys the 
su(2) <8> sw(2)-algebra. The original form of the Hamiltonian (12- lj) can be expressed in terms 
of fermion operators shown in the form (1-2-1). The explicit expression of the Hamiltonian 
(12-11) in boson operators is given in the relation (I-3-3b). 

For the present boson system, the following four hermitian operators are mutually com- 
mutable: 

L = (h/2)(a* + a+ + b* + 6+) , (2-3a) 
M = {h/2)(a* + a + + a*_a_) , (2-3b) 
f = (h/2)(-a* + a + + b*_b_ + 1) , (2-3c) 
K = ha* + a + . (2-3d) 

The form (I2-3P is shown in the relation (1-3-4) and it should be noted that L, M and T 
commute with the Hamiltonian (12-ip . With the use of the operators (12-31) . together with the 
expression (!2-2[) . the Hamiltonian (12-11) can be rewritten as 

H = % + H 1) (2-4) 



Hn 



Hi 



(L + M - (T - h/2^j + 4GTM 



+2 



e-G [L + M - [T- h/2 



K + 2GK 2 



-G-h 2 



a*,b*_b + a_ + a*_bVb-Ci A 



(2-5a) 
(2-5b) 



In (I), noting that L, M and T are constants of motion, we treated the Hamiltonian f)2-4p 
quantum- and classical-mechanically. In the case of the classical treatment in (I), the sw(2)® 
su(l, l)-coherent state played a central role. In this paper, we treat the cases of the su(l, 1)® 
su(l, 1)- and the su{2) ® sw(2)-coherent states on an equal footing with the su{2) (g> su(l, 1)- 
coherent state. 

As was mentioned in (I), the two-level pairing model as a many-fermion system is char- 
acterized by three quantities. They are the numbers of the single-particle states in the levels 
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a = + and a = — , namely, 2i? + and 2i?_, respectively, and the total fermion number N. 
The original fermion system is related to the present boson system through 

(H/2)a+ = L, (h/2)N = M, {h/2){Q- — N/2 + 1) = T . (2-6) 

Here, the g-numbers L, M and T are replaced with the c-numbers L, M and T, respectively. 
Through the relation ( 12-61) . the subspace, in which the eigenvalues of L, M and T are L, M 
and T, respectively, corresponds to the original fermion space specified by f2 + , Q- and N. 
The above can be found in the relation (1-6-3) . In this paper, we treat the case T > h/2. 
For the numerical results, we discuss the case 

Q + = Q_ (= O) , N = 2Q . (2-7) 

The case (12-71) corresponds to 

L = M = (h/2)Q , T = h/2 . (2-8) 

The case (12-71) tells that the degeneracies of the two levels are the same and the total fermions 
occupy the lower level completely (a closed shell), if the interaction is switched off. 

§3. Re-formation in terms of one kind of boson operator and its classical 

counterpart 

The system presented in §2 can be expressed in terms of four kinds of bosons and it 
contains three constants of motion. Therefore, the present system can be described in terms 
of one kind of degree of freedom. With the aid of the MYT mapping method, in (I), we 
transcribed the system in the space spanned by one kind of boson (c, c*). The Hamiltonian 
(12 -4p is transcribed in the form 

H = H + Hi, (3-1) 



H = - [e (L + M - (T - h/2)) + AGTM] 

+2 [e - G (L + M - (T - h/2))] K + 2GK 2 
H 1 = H 1 {c^) = H 1 {c i) )=H 1 {c) , 

HJc) = -G 



(3-2a) 



Vhc* ■ \ 2T + K\ 2M - K\ 2L - K 



+ \ 2L- K\ 2M - K\ 2T + K-Vh 
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Mc ) = -G 



^2L + h-k ■ Vhc* ■ \Jit + k\jiM -k 



H x {c) = -G 



+^2M - K\jlT + K -Vhc- \j2L + h-K 
2T - h + K ■ Vhc* ■ \]2M - K\]2L- k 



+ yj2L- K\j2M - K -Vhc- \j2T-h + K 



K = hc*c . 



(3-2b) 
(3-3) 



The notations c_, c° and c will be interpreted later. The relation [^/hc* , K] = —h - Vhc* 
and [Vhc , K] = h - Vhc shows that i?i(c_) = Hi(c°) = H(c). The above can be seen in the 
relations (1-4-11) and (1-4-12) with the derivation. 

We know that, with the aid of the appropriately chosen boson coherent state, we can 
derive the classical counterpart of the original quantal system. For the present system, we 
used in (I), the coherent state |c_)) shown in the form (1-5-3) and its rewritten form (1-5-6). 
The state |c_)) can be, further, rewritten as 



|c_)) = N c exp 



U. 



+ 



x exp 



Y± 

A- 



■a*_b*_ ] exp 



J_W- 

vnTT + 

a*_b + ) exp ( -^B + b* + ) |0» , 



V + 



(3-4) 



Here, A^ c denotes the normalization constant and V+, A- and B + are complex pa- 

rameters. We can see that the state |c_)) is generated by successive operation of a*_b + 
and a* + b*_. The operators a*_\ b + and a* + b*__ are the raising operators of the su(2)- and the 
su(l, l)-algebra, respectively. From the above reason, in (I), we called the state |c_)) the 
su{2) <S> su(l, l)-coherent state. Through the process presented in (I), the expectation value 
of the Hamiltonian if, which is given in the relation (1-3-3), was calculated. The result is 
shown in the relation (1-4-15). In notations slightly different from those in (I), we have the 
following form: 



if(c_) = ({c-\H\c-)) = H (c-) + ifi(c. 



(3-5) 



H (c-) — — [e(L + M — (T — h/2)) + 4GTM] 

+2 [e - G (L + M - (T - h/2))] K + 2GK 2 , 
H^c-) = -2GVK(2T + K)(2M - K)(2L - K) cos^ - 



(3-6a) 
(3-6b) 
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Here, (ip, K) denotes the angle-action variable. 

In (I), we showed various forms of coherent states. One of them was presented in the 
form (1-7-20), which is denoted as |c )). The state |c )) can be rewritten as 

W = ^-«p(^fc)exp(^fc) 

x exp ex P |0)) ' 

U a = Vl + IKI 2 - (3-7) 

Here, of course, N c denotes a normalization constant and V a and W a (<r = ±) are complex 
parameters. We can see that |c )) is generated by the operators a* + b*_ and a*_b* +) which are 
the raising operators of two independent su(l, l)-algebras. Then, we call the state (13 -TP 
the sw(l, 1) <g> su(l, l)-coherent state. Under the same idea as that in the case of |c_)), the 
expectation value of H for |c )) is given in the form 

H(c°) = ((c°\H\c )) = H (c°) + Hxic ) , (3-8) 



H (c°)=H (c_) , (3-9a) 
H^c ) = -2G^K(2T + K)(2M -K)(2L + H- K) cosip . (3-9b) 



The third is the following one: 



iV c exp ( — -alS_ J exp ( —=B_b* 



x exp (^2-6+) ex P (-^ B + b *+) l°» • ( 3 ' 10 ) 

Here, iV c denotes the normalization constant and v4 CT and _B CT (cr = ±) are complex param- 
eters. Clearly, |c)) is generated by a* + b^ and a*_b + , which are the raising operators of two 
independent sw(2)-algebras. Therefore, we call it the su{2) ® sw(2)-coherent state. We did 
not contact with |c)) explicitly in (I). The state |c)) can be rewritten as 



|c)) = iV c exp ( ^±a* + b + J exp (^=&+ 



x exp ( ^a*J_ j exp f ^6* ) |0)) 



iV c exp f (^+0+ + 5+6; + A_al + |0)) . (3-11) 
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The form (13-111) tells that |c)) corresponds to the BCS state and it is identical to a kind of 
the Glauber coherent state. Then, the expectation value of H for |c}) is calculated in the 
form 

H(c) = ((c\H\c)) = ff (c) + H^c) , (3-12) 
H (c) = F (c_) , (3-13a) 



#i(c) = -2G^K{2T -h + K){2M - K)(2L - K) cos ip . (3-13b) 

In the above, we showed that, for one quantal system, three classical counterparts were 
derived. It may be natural, because we used three different coherent states. Our final 
problem of this section is related to the re- quantization. For this task, it may be convenient 
to introduce new canonical variable in boson-type, (c, c*), through 

Vhc = VKe'^ , Vhc* = y/Kj* . (K = hc*c) (3-14) 

The re- quantization may be performed by the replacement 

Vhc — ► Vhc , Vhc* — > Vfre* , K — >K . (3-15) 

Under the replacement (13-151) . H Q (cJ) = H (c°) = H (c) is requantized to H shown in the 
relation (l3-2al) . Concerning Hi(c-), Hi(c°) and Hi(c), the following relation is interesting: 



VhcW2T + KV2M - KV2L - K — > Vhc*\/2T + K\J2M - K\J 2L - K 

for |c_» (3-16a) 
y/hc*V2T + KV2M - K^2L + h-K = ^2L + h- kVHcW2T + K^2M - K 



2L + h-kVhc*\/2T + k\/2M-k = Vhc*\/2T + k\]2M - k\]2L- K 

for |c )) (3-16b) 
Vhc*V2T -h + KV2M - KV2L - K = \J2T -h + KVhc*V2M - K^2L - K 



2T-h + kVhc*y 2M - ky 2L - k = Vhc*y2T + ky 2M - ky 2L-k 

for |c)) (3- 16c) 

If we note the relation 2K cos?/' = \/K(yhc* + Vhc), the relation (13-161) leads to 

H x {cJ) — . F x (c_) , ^(c°) — . ^(c ) , H x {c) — . H x {c) . (3-17) 



As was shown in the relation (13-2bj) . we have H\(cJ) = Hi(c°) = H\{c) = H\. In the above, 



the meanings of the symbols c_, c° and c may be clear. From the above argument, we can 
understand that three different classical systems become to single quantal system, depending 
on the ordering of the variables. Therefore, the above three are on equal footing, and then, 
it may be interesting to investigate what results they give for each classical and quantal case. 
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§4. Classical treatment 



In §7 of (I), we sketched an idea how to treat the Hamiltonian obtained in classical and 
quantum framework. In this section, we apply this idea to the three classical cases discussed 
in §3. We treat the concrete case (12-81) which is equivalent to the condition (12-71) in the 
original fermion space. Further, in order to simplify various relations, we adopt h = 1 and 
e = 2. Then, the functions F(K) and f(K) introduced in the form (1-7-1) is given in the 
following form: 

H = F{K) - f(K) + 2f(K) (sin^/2) 2 . (4-1) 
F(K)=F + F 1 (K) , (4-2) 



F = -{2 + G)Q , 

F X {K) = 2(2 - GQ)K + 2GK 2 



f(K) = 2G V / K(K + 5)({2-K)(f2 + A-K) (= f(K; 5, A)) 
Here, (5, A) denote parameters for discriminating the three cases: 



(4-3a) 
(4-3b) 

(4-4) 



5 = 


1 , A 


= 0) 


for 


#(c_ 


5 = 


1, A 


= 1) 


for 


H(c° 


5 = 


, A 


= 0) 


for 


H(c) 



(4-5) 

Our picture is based on small amplitude oscillation around the energy minimum point. 
If the energy minimum point can be found in the region < K < Q for G ^ 0, the form 
flUU gives 

ip = . 



Then, our problem is reduced to finding a solution of the relation 

F\K) - f(K) = , i.e., F[{K) - f{K) = 

Here, F[(K) and f'(K) are given as 

F[{K) = 2(2 - GO) + AGK , 
2K + 5 



(4-6) 



(4-7) 



f'(K) = G 



-y/K(K + 6) 



y/(f2 - K){Q + A- K) 



2Q + A-2K 



y/(n-K)(fi + A-K) 



(4-8) 
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Fi(K) 




(a) (b) 



Fig. 1. The behavior of F\(K) is depicted in the case (a) GO, < 2 (i? = 5 and G = 0.25) and (b) 
GQ > 2 (J7 = 5 and G = 0.75). Here, ^'(0) = 2(2 - GQ). 



f (K) 




Fig. 2. The behavior of f(K) is depicted in the cases of (a) 5 = 1, A = (case (i)), (b) 5 = 
1, A = 1 (case (h)) and (a) 5 = 0, Z\ = (case (hi)). In (a) and (b), /'(0) goes to infinity, 
/'(0) -> oo. On the other hand, in (c), f'(0) = 2GQ. 



Figs.l and 2 show that F{(0) is finite and /'(0) — >■ +oo in the case 5 = 1, and (Fi(K)-f(K)) 
is decreasing near K = and then, afterward, changes to increasing function until K — Q. 
Therefore, we can find the solution of the relation (14-71) in the region < K < Q. On the 
other hand, in the case (6 = 0, A = 0), if ^'(0) - f (0) > 0, i.e., GQ < 1, {F X {K) - f(K)) 
is an increasing function in the region < K < fi and we cannot find the solution of the 
relation (14-61) . In this case, K = gives the energy minimum point and the condition (14-60 
is not necessary. If F{(0) - f'(0) < 0, i.e., GQ > 1, (F^K) - f(K)) is decreasing near K = 
and then, afterward, changes to increasing function. Therefore, we can find the solution of 
the relation f)4-6p . The above consideration gives the energy of the energy minimum point 
in the form 

E C (K ; 5, A) = F + F.iKo) - f(K ; 5, A) . (4-9) 

Here, for the cases 5 = 1 and 5 = (GQ > 1), K is given as the solution of the relation 
(EU): 

F[(K )-f\K ;5,A) = 0. (4-10) 
For the case 5 = (GQ < 1), K is given as 

K = . (4-11) 
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The frequency 10 which characterizes the small amplitude oscillation around the energy 
minimum point is given in the form for the cases 5=1 and 5 = (GQ > 1) 



co = [(F"(K )-f"(K ))f(K )] 1/2 . 



(4-12) 



This form can be found in the relation (1-7-5). Of course, the frequency (14-121) depends on 
K , 5 and A. Then, denoting as uj(K ; 5, A) and using the relation (14-101) . we have 

1/2 



oj(K ;5,A) 



F[(K Q ) 2 + F?(K )f(K ; 5, A) - (1/2) (/ 2 )" (K ; 5, A) 



(4-13) 



Here, (f 2 )"(K; 5, A) denotes second derivative of f(K;5, A) 2 for K. The frequency in the 
case 5 = (GQ < 1) is also obtained. The method is discussed in the Appendix B. Thus, 
the Hamiltonian H is expressed as 



H = E C (K ; 5, A) + co(K ; 5, A)d*d . 

Here, (d, d*) denotes boson-type canonical variable. 

The above general argument presents the following concrete expressions: 
(i) For |c_» (5 = 1, A = 0): 



(4-14) 



1 0) 2l \K + 1 4y/K (K + 1) 



w(c_) 



-2Q + AK -G 



Q + 2 v / ^o (v^o + VKq + l) (Q - K ] 
(2 - GQ) 2 + 4G(2 - GQ)K 

G 2 (2K 2 - 3(2Q - 1)K + Q(Q - 2) - 2^K Q (K Q + 1)(Q - K ] 



1/2 



(4-15) 



(ii) for |c )) (5 = 1, A = l] 



G 



-1 



(Q - 2K ) - - I 1 



+ 



(^0 + 1)^-^0 + 1) 

(fi-2ifo)(fi + l) 



^c(c°) 



iy/K (K + - # )(^ -K + l) 
-2Q + 4K 



-G 



Q + 2v/^ (^ - ATo) (y/K ({2 - tf ) + v^o + 1)(Q-K + f 



cu(c°) = 2 



(2 - Gi?) 2 + 4G(2 - GJ2)K 
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-30 



0.1 0.2 



0.3 



0.4 

G 



0.5 0.6 0.7 



0.8 



Fig. 3. The classical energies are depicted as a function of G in the cases of (i) 5 = 1 and A = 
0E? c (c_)), (ii) 8 = 1 and A = 1 {E c (c )) and (iii) 5 = and A = (£ c (c)), where ^ = 5. Here, 
-Eexact shows the exact eigenvalue of the Hamiltonian. 



-G 2 (2K 2 - QQK + (Q 2 - Q - 1) 



1/2 



(4-16) 



(iii) for |c)) (5 = 0, Z\ = 0, GQ < 1): 



K = , 

£ c (c) = -2Q - GQ 
uj{c) = 4Vl - GQ , 



(4- 17a) 



(iii)' for |c)) (5 = 0, A = 0, GM? > 1): 



G" 1 = J? - 2^0 , 

£ c ( c ) = cr 1 - GQ{Q + 1) 

w(c) = 2y/(GQ) 2 - 1 . 



(4- 17b) 



Here, each i? c is identical to the energy in the classical description of the two-level pairing 
model. In this paper, we have presented three different classical counterparts for the original 
system, which lead to the corresponding energies, _E c (c_), E c (c°) and E c (c), respectively. 
Figure [3] shows them together with exact eigenvalue E exact . Then E c (c) reveals the energy 
obtained by using the su(2) g) s-u(2)-coherent state, that is, the Glauber coherent state, 
which corresponds to the energy obtained by the Hartree-Fock-Bogoliubov calculation in 
the original fermion system. It seems that the use of the su{2) <g) su(l, Incoherent state, 
|c_)), and the su(l, 1) ® su(l, l)-coherent state, |c )), comparatively give good results in 
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comparison with the exact energy eigenvalue in the region of small force strength G. On the 
other hand, the su(2) ® s-u(2)-coherent state |c)) gives rather good result in the region of 
large force strength G in this stage. In the next section, we reinvestigate the ground state 
energy and the oscillation frequency around the static configuration in terms of the quantal 
description based on these three classical counterparts. 

§5. Quantal treatment including quantum fluctuations 

In order to compare the energies derived in the above procedure with the exact ground- 
state energy in quantal description, we introduce the energy E q with quantum correction 8 
as 



E q = E c + £ 
E = \ 



f(K ) g'(K 0) 



K o g{Ko) 

where g(K) = g(K] 5, A) is defined through f(K) in (j4-4p as 



(5-1) 



f(K; 5, A) = 2G0f(FTT)(« - K)J 



(5-2a) 



g{K;5 = 


1,A = 


0) 


= 1, 


g(K;S = 


1,A = 


1) 


= Q-K + l , 


g(K;6 = 


0,A = 


0) 


= K . 



(5-2b) 

In a form similar to the Hamiltonian (I4-I4p , we have the quantized Hamiltonian in the form 

H = E q (K ; 5, A) + uj(K ; 5, A)d*d . (5-3) 

Here, it should be noted that the frequency uo, which is given in the relation f)4-13j) . plays 
the role of the excitation energy in the present case. The derivation of f!5-3p with (15-11) and 
the frequency ui is given in appendix B. 

Figure H] shows the behavior of E in the cases (i), (ii) and (iii), compared with the exact 
ground-state energy E^^, as a function of the strength of interaction G. The behavior is 
almost the same in each case except for case (iii) near G = 0.2. In the case (iii), a dip appears 
at G = 0.2, which corresponds to the phase transition point in the Hartree-Fock-Bogoliubov 
calculation, namely, GQ = 1 with Q = 5. The su(2) <8> s-u(2)-coherent state |c)) in (13-111) in 
case (iii) is identical with the Glauber coherent state. This state corresponds to the Hartree- 
Fock and/or the Hartree-Fock-Bogoliubov state in the original fermion space. This behavior 
near G = 0.2 in case (iii) is realized in the quasi-particle random phase approximation 
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_30 I i i i i i i i I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
G 

Fig. 4. The energies are depicted as a function of G in the cases of (i) 6 = 1 and A = (-EVn), (h) 
6 = 1 and A = 1 (-Enn) and (iii) 5 = and Z\ = (-Bpy)), where Q = 5. Here, E exa _ ct shows the 
exact eigenvalue of the Hamiltonian. 
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2 - V~"" /" 

1 " 

o I i i i i i i i 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



G 

Fig. 5. The frequencies are depicted as a function of G in the cases of (i) 5 = 1 and A = 

(ii) 6 = 1 and A = 1 (wpj) and (iii) 5 = and Z\ = (w^), where Q = 5. Here, <^ cxa ct shows 
the exact result. 

(QRPA). n) Thus, it is seen that the energy correction 8 gives the ground state correlation 
in the random phase approximation. On the other hand, in the cases (i) and (ii), the dip 
structure does not appear. 

Figure [5] shows the behavior of the frequency uj in the cases (i), (ii) and (iii), compared 
with the exact frequency cj exact , which is exactly calculated for the excitation energy from the 
ground state. In the su{2) (g> sw(2)-coherent state |c)), the dip structure appear at G = 0.2, 
which is the same behavior as the QRPA calculation. However, by using the su{2)®su(l, In- 
coherent state |c_)) in case (i) and the su(l, 1) <g>su(l, Incoherent state |c )) in case (ii), this 
fictitious dip structure does not appear and the derived frequencies well reproduce the exact 
frequency. Especially, the result in the case (i) shows same behavior as the exact frequency 
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in the region below the appearance of the dip in case (iii). However, the case (iii) based on 
the su{2) (g> sw(2)-coherent state approaches the exact frequency in the region of the large 
G. 

§6. Concluding remarks 

We investigate the ground-state energy and the frequency of the small amplitude oscilla- 
tion around the ground state in the two-level pairing model governed by the su(2) <8> su{2)- 
algebra with the viewpoint of using the various states based on the su(2)- and the su(l, 1)- 
algebras in the boson representations. Various classical descriptions are possible through the 
various states mentioned above, while the original fermionic quantum theory is unique. 

The su(2) ®s«(2)-coherent state in the boson representation is identical with the Glauber 
coherent state, which corresponds to the BCS state (Hartree-Fock-Bogoliubov state) in the 
QRPA calculation. Then, the fictitious sharp phase transition point appears. However, 
taking into account the possibility of the various classical descriptions based on the various 
states, the fictitious sharp phase transition can be avoided. Thus, a reasonable approximation 
can be obtained by using other coherent states such as the su{2) ® su(l, 1)- and su(l, 1) <8> 
su(l, Incoherent states. 

Further numerical results will be reported in a subsequent paper, where g(K) introduced 
in the relation fl5-2al) is adopted in a form different from those shown in the relation f!5-2bj) . 
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Appendix A 

Proof of the result ^TTTa ) 

In this appendix, we show a proof of the result (I4-I7al) . For this aim, let us start with 
the classical Hamiltonian (13- 12j) . which was obtained under the Glauber coherent state, and 
treat the case in which the effect of K is negligibly small compared with hQ. Under this 
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condition, by picking up linear terms for K, the Hamiltonian (13-121) is approximated to 



H = -h(e + hGQ) + 2(e - hGQ)K - 2hGQK cosip . (A-l) 

Of course, we consider the case L = M = hfi/2 and T = h/2. Since (ip,K) is a set of 
canonical variables, we have the following Hamilton equation of motion: 

Tp =0K= 2 (e- hGQ) - 2hGHcosip , (A- 2a) 

K=-— = -2KGQK sini) . (A-2b) 
dip 

We can solve Eq.( 1A-2l) directly. But, it may be much more transparent to solve it by 
introducing the quantities 

T x = Kcosip , T y = K sin ip , T z = K . (A-3) 

These are constrained by the condition 

T |_( T 2 + r y 2 )=0. (A-4) 

We describe the present system in terms of (T x ,T y ,T z ). With the use of the relation (1A-3I) . 
the time-derivative (T x ,T y ,T z ) is given in the form 

f x = -2(e - hGQ)T y , (A-5a) 

f y = 2(e - HGQ)T X - 2hGf2T z , (A-5b) 

f z = -2hGQT y . (A-5c) 



The relation (1A-5P gives us 



T y = -{2ef ( 1 - ^) T y . (A-6) 



Since T y should be finite, the solution of Eq. (IA-6j) is meaningful in the case 

2hGQ 



> . (A-7) 



Then, we define the frequency uj in the form 



2hGQ /A 
u = 2e\l . (A-8) 
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It may be clear that T y shows a harmonic oscillator type behavior and the relation ( 1A-4I) 
and (1A-5j) give 



T„ = A 



2U-hGQ) , , 2/lGJ? 
- cos(cjt + a) + 



Ty = A sm(ut + a) , 
T r = A 



IhGQ , . 2(e-frGJ?) 
cos(c<ji + a) + - 



a; 



(A-9a) 
(A-9b) 
(A-9c) 



Here, A and a denote constants to be determined by the initial condition. The Hamiltonian 



(1ATI) reduces to 



H = -H(e + HG)Q + ujA . 



(A-10) 



If we put h = 1 and e = 2, the above result obtains the form (I4-I7al) . Later, the meaning of 
D will be discussed. 

If we disregard the starting condition K UQ, the present result can be accepted in 
the region 2HGQ /e < 1 and any value of A. If we keep in mind the starting condition, some 
discussions are necessary. Since — 1 < cos# < 1 and T z = K, the form ( 1A-9cl) gives us 



(A-ll) 



Then, combining K <C fiQ with the relation (lATip . we have A(2e/u) <ti hf2, i.e., 



A<g:hf2 — . 

2c 



(A-12) 



If becomes smaller, the value of A may also become smaller. 

It may be interesting to see that (T x ,T y ,T z ) obeys the su(l, l)-algebra classically. This 
can be seen in the relation 



T + , T_ }p = 2iT n 



[T ,T±)p = TiT± . 



(A-13) 



Here, [ , }p denotes the Poisson bracket and (T± i0 ) is defined as 



T+ = T X ± iT n 



(A-14) 



The relation (1A-4I) is rewritten as 



T 2 - T + T_ = . 



(A-15) 
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Appendix B 

The derivation of quantal Hamiltonian with the disguised form 



In this appendix, we give the quantal Hamiltonian describing the small oscillation around 
the energy minimum state in a slightly disguised form for the original Hamiltonian (13- II) . 
A part of the Hamiltonian, Hi, is recast into 



-G 



-G 



Vhc* ■ \ 2T + K\ 2M - K\ 2L-K + h.c. 



(A) ( Vhc* - V k) A/2T + kJ 2M - A A/ 2L — K 



g(K+h) 



+ VK\ 2T + K\ 2M- K\ 2L - A 



9(K) 
g(K + h) 



-h.c. 



;b-i) 



where A = hc*c and g(k) is an arbitrary function with respect to A. Here, h.c. represents 
the Hermite conjugate. Thus, the Hamiltonian (13- ip can be expressed as 



H = F + Fx(K)-f(K) 

-- 2G\J g{K)(Vhc* - \fk)\j2T + K\J 2M — K\J 2L — A = 

21 \Jg(K + h) 

F = -[ e -G(L + M - (T - h/2)) + AGTM] , 
Ax (A) = 2 [e — G(L + M — (T — H/2))} A + 2GA 2 , 

9(K) 



+ h.c. 



/(A) = 2GV Ay 2T + Ky2M- Ky2L - A 



s(A+ft) 



;b-2) 



Here, let us focus on the last bracket in the Hamiltonian H. From (I- 7-11), we know that 
Vhc* - \fk = iVh(-y* + j)/2 - ((V^(7* + 7)) 2 - 2fr)/8 v / Ao. Then, substituting the above 
relation into the last term and expanding the terms up to the second order of (7,7*), we 
obtain 

H = F + Fx(K ) - /(A ) + (F{(K ) - f\K ))6K + ^(Af (A ) - f'(A ))(5A) 2 



7=^+7)) + l(f(k 



/(A ) 9'(K ) 
0) 7? tttt/I^o, 



An 



g(K 



(B-3) 



Here, we have introduced 5A = A — A = v^A) • iVhij* — 7) + ^7*7 in (1-7-8) and used 
the commutation relation [ <5A , iy/h^* + / y)/2 ] = H^Kq. Then, Ao is determined by 
eliminating the linear term of <5A as F[(Ko) — /'(An) = in Eq. (14 -7p . Further, the above 
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mentioned commutation relation can be expressed as [ 5K , v^(7* + 7)/2a/^o ] = —ih. 
Thus, the above Hamiltonian (1B-3j) is easily diagonalized by means of newly introduced 



operator (d, d*), which is written as d — u\/h{ A f* +^f) + v6K with u-v factor satisfying 
u 2 — \v | 2 = 1. We finally obtain the Hamiltonian as 



H = F + F^Kq) - f(K ) + hwd*d 



CO 



V(F{'(K ) - fMpoj • (B-5) 



In this paper, since we discuss the case L = M = {h/2)fi and T = h/2 in Eq. (12-81) and the 
parameter e = 2 in the natural unit H — 1, we obtain (14- 3 p and (15- 21) . 
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